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Abstract 

We apply the dinamically driven renormalization group to study 
the critical properties of the Domany-Kinzel probabilistic cellular au- 
tomaton. To preserve the absorbing state clusters with at least one 
site occupied are renormalized into a one occupied site and clusters 
with all sites empty are renormalized into a one empty site. We have 
obtained the phase diagram as well as the critical exponent related to 
the spatial correlation length. 
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1 Introduction 



Among the variety of irreversible models that describe a nonequilibrium 
phase transition from an active state to an absorbing state, the Domany- 
Kinzel probabilistic cellular automaton [|IJ together with the contact process 
are perhaps the simplest ones. According to Janssen || and Grassberger 
0], any one component model displaying a transition from an absorbing to 
an active state belongs in the same universality class as direct percolation. 
In fact, the Domany-Kinzel cellular automaton is equivalent to direct perco- 
lation. 

The Domany-Kinzel probabilistic cellular automaton is defined as follows. 
In a one dimensional lattice each site can be in one of two states, either 
empty or occupied. At each time step the sites are updated simultaneously 
and independently. The state of a given site depends only on the states of 
two neighboring sites in the previous time step. The rules are: 1) if the two 
sites are both empty, the given site will be empty, 2) if one site is occupied 
and the other is empty, the given site will be occupied with probability pi, 
and 3) if the two sites are both occupied, the given site will be occupied with 
probability p2- Due to the rule number 1, the state with all sites empty will 
be an absorbing state. 

In the stationary state the Domany-Kinzel automaton exhibits two states: 
the absorbing state with all sites empty and an active state with a nonzero 
density of occupied sites. For small values of pi the system is in the absorbing 
state. When the parameter p\ is increased the absorbing state becomes 
unstable, for sufficiently large values of pi, given rise to an active state. 
The phase transition from the absorbing to the active state is a continuous 
transition and belongs to the same universality class as direct percolation. 

Renormalization group (RG) techniques in real space have been used 
successfully in studying the critical behavior of equilibrium statistical mod- 
els since the work of Niemeyer and van Leeuwen ||. In these techniques 
small clusters of sites are renormalized into just one site according to certain 
rules. For the case of the Ising model with the up-down symmetry, the most 
important rule is the majority rule which preserves the up-down symmetry 
of the model at each step of the renormalization. Any RG scheme should 
be set up in a way to preserve the symmetries and the main properties of 
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the model studied. In the present case the most important property is the 
existence of an absorbing state which should then be preserved by the RG 
trasformation. To this end, we have constructed a real space RG such that 
a cluster with at least one site occupied is renormalized into a one occupied 
site. Only clusters with all sites empty are renormalized into a one empty 
site. 

For systems in equilibrium, that is, for systems defined by a Hamilto- 
nian the appropriate RG space is the space of the parameters, or coupling 
constants, that define the Hamiltonian. A sequence of RG transformations 
correspond to a trajectory in this space. Systems far from equilibrium, that 
is, systems that lack detailed balance, on the other hand, are not defined 
by Hamiltonians but by the dynamic rules or equivalently by the transition 
probabilities that govern its time evolution. In this case the appropriate RG 
space is the space of parameters that define the transition probabilities. 

Here we use a real space RG || |§ that renormalizes the transition 
probabilities so that one obtains a RG transformation in the space spanned 
by the parameters that define the transition probabilities. In the case of the 
Domany-Kinzel cellular automaton studied here this space will be the space 
defined by the parameters pi, and p%- The scheme we use is an application 
of the dynamically driven renormalization group introduced by Pietronero et 
al. and Vespignani et al. |TD[ [[□]] to study systems in critical states. A 
study of the Domany-Kinzel by means of a phenomenological RG has been 
done by Bagnoli et al. JT2|. 

The RG transformation was accomplished by using cells of size 6 = 2 and 
by using a series of approximations to the stationary state. As we increase 
the order of the approximation better values were obtained for the universal 
critical exponent as well as for the nonuniversal critical quantities. Our best 
value for the critical exponent uj_ related to spatial correlation length is 
u± = 1.013 which should be compared to u± = 1.100 ± 0.005 obtained from 
transfer matrix methods [|IJ. The critical value of pi along the line P2 = 
is p c = 0.792 which is close to p c = 0.799 ± 0.002 obtained from simulations 



2 RG transformation 

Consider a one dimensional lattice of N sites in which each site can be in 
one of two states, either empty (dj = 0) or occupied (cr^ = 1). At each 
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time step the sites are updated simultaneously, and independently, so that 
the transition probability VK(cr|cr') from state a' = (a[,o J 2 ,...,a' N ) to state 
a = (di, cr 2 , cr N ) is given by 

W(a\a')=Uw(a l \a' l ,a' l+1 ), (1) 

i 

where 

w(l|00) = 0, io(l|10) = w(l|01) =px, w(l|ll)=p2- (2) 

These rules completely defines the Domany-Kinzel probabilistic cellular au- 
tomaton since 

w(0\a>,a' l+1 )+w(l\a>,a' l+1 ) = l. (3) 

We remark that the expression implies that, given the state cr', the 
variables oi are independent random variables which in turn is equivalent to 
say that each site should be updated independently of the others, as desired. 
We call this property of the cellular automaton one-site independence. 

In the stationary state the following condition 

P(a) = J2W"(a\a')P(a') (4) 

cr' 

holds for any value of £, where P(o~) is the stationary probability. If we use 
the notation Pe(a, cr') for the probability of having state cr' at a given time 
and state cr at a £ time steps later then 

P e (a,a')=W e (a\a')P(a'). (5) 

Similarly, 

P(r,r') = W(t\t')P(t') (6) 

is the analog equation for the renormalized system where W(t|t') is the (one- 
step) renormalized transition probability we are looking for and r = (ti,T2, 
...,Tjv')) with N' = N/b, denotes a renormalized state. 

Let TZ(r\a) be a conditional probability of the renormalized state r given 
the state a, which has the following properties 

^(r|a)>0, 52K(t\<t) = 1. (7) 
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The RG transformations are obtained by demanding that || 

Hr,r) -EE^k)^V0^(^O (8) 

<r cr' 

from which follows 

P( t ')=£^(tV)^V)- (9) 

This last expression is the one used in real space RG for equilibrium systems. 
Taking the ratio of the last two equations we obtain |TT] 





a)W l { 


a 


a')n(r' 


o')P{o>) 




a 


>)P{<j<) 





W(t\t') = v ' 'V'^ ' ' V , (10) 



which is the desired expression for the renormalized transition probability. 
Due to the presence of the stationary probability on the right hand side of 



(|T0|), this is not a transformation that involves only the transition probability. 



However, if we use an approximation for the steady state P(o~), such as 



a mean field approximation, then equation (10) provides a well defined RG 
transformation W — *> W from the old W to the new W transition probability. 
It is interesting to write equation fllCf ) in the form 

w(t\t') = j2Y, n ( T \ a ) wi ( a \ a ') s ( a '\ T ') ( n ) 

a cr 1 

where 

n(T'\e')P(a') 
SHT)= £„.R(rV)P(ar <12) 
From this last expression it is clear that S(a'\r') is normalized with respect 
to the variables a' and that it can be interpreted as a conditional probability 
of a' given r'. In this way it becomes quite obvious that the renormalized 
transition probability W(t\t') is properly normalized with respect to the 
variables r as it should. One important property of S follows from a specific 
property of TZ. We will see in the next section that lZ(r\cr) is either zero or 
one from which follows that 1Z(t' \<j)1Z(t\<j) = 8{t' ,t)1Z(t\ct). Using this last 
equality we can check that 

Y j n{T\a)S{a\r') = 5{r\T), (13) 

17 

independently of -P(cr). 
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We remark that we should not expect, in general, that, given the state 
r', the variables Tj be independent. In other words, the stochastic process 
defined by W(t\t'), given by equation (|I0D, will not have the desired one-site 
independence. This property, however, can be accomplished in an approxi- 
mate way as follows. Let us define the one-site marginal transition probability 
w{ji\T') by 

w(t 1 \t>) = J2'W(t\t>), (14) 

T 

where the sum is over all Tj except Tj. Next we construct a transition proba- 
bility Wr(t\t') as being the product of these one-site transition probabilities, 
that is, we define 

W r {t\t') = X[w{t 1 \t'). (15) 

Of course Wr{t\t') does not coincide, in general, with W(t\t'). 

The final step in the renormalization procedure is to say that the renor- 
malized cellular automaton will be defined by the transition probability 
Wr(t\t') which in turn has the form of equation ([!]). The last step is an 
approximation which we call the one-site independence approximation. In 
this manner we obtain a renormalized probabilistic cellular automaton, with 
the property of one-site independence, whose one-site transition probability 
w(Ti\r') is obtained from the original one w(ai\a') through equations (|i~4"D, 
(|T0|), and (0). Thus we understand the present renormalization as the trans- 
formation w (<7i\cr') — > w{Ti\r'). 

3 Renormalization algorithm 

We use a renormalization 1Z(t\<j) in the form of a product which for the case 
6 = 2 reads 

N/2 

TZ(r\a) = Y[ R{Tk\°2k,V2k+x), (16) 
fe=i 

where 

i2(l|00) = 0, #(1110) = #(1|01) = i?(l|H) = 1, (17) 

and 

i?(0|00) = 1, i?(0|10) = i?(0|01) = i?(0|ll) = 0. (18) 

The rule _R(0|00) = 1 ensures that the state with all sites empty will be again 
an absorbing state after the renormalization transformation. 
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To calculate the renormalized transition probability w(ti\t2,t 3 ) we con- 
sider a cluster of nine sites of the original space-time lattice: two sites (labeled 
1 and 2) at a time £, three (labeled 3, 4 and 5) at time i — and four (labeled 
6, 7, 8 and 9) at time 1 — 1. The renormalization of the sites are as follows: 
(0-1,02) -> ri, (0-6,0-7) -» r 2 , and (08,09) -> r 3 . Using equations (|l|), (|12D, 
and (0), the renormalized transition probability w{ti\t2,t 3 ) are calculated 
by 

w{ti |r 2 , ^)=53S53SSS -ftfaki* ^2)^2(^1, 02K, 07, o 8 , (J9) X 

CTi (72 0-5 0-7 erg erg 

x p(or 6 ,cr7,o- 8 ,cr 9 |T2,T3), (19) 

where 

102(01, 2 |cr 6 , 07, 8 , Og) = 

= mi E W ( ff ll ff 3. 04)^(02104, 5 )w(o 3 |o 6 , 7 )w(0 4 |07, 8 )w(o 5 |o 8 , Og), 
0-3 0-4 05 

(20) 

and 

p(cr 6 , 07, o 8 , o- 9 |r 2 , r 3 ) 
with 

P(r 2 ,r 3 ) =^2^2^2^2R(T 2 \ae,a 7 )R(T 3 \a 8 ,a 9 )P(a 6 ,a 7 ,a 8 ,a 9 ). (22) 

""6 C"7 erg erg 

To actually use the above equations we have to calculate in each step 
of the renormalization the stationary probability distribution -P(a) related 
to the old transition probability w{ai\a' il o- +1 ) which satisfies the balance 
equation 

P(a)=Y,W(a\a')P(a>), (23) 

cr' 

where ^(ojo') is given by ([!]). Of course, this equation cannot be solved 
exactly and we should seek for approximate solutions. 

The simplest approximate solution of (|23|) is obtained by writing the 
following equation for the one-site probability 

P(^) = J2J2w(a 1 \a' 1 ,a' 2 )P(a[,a' 2 ), (24) 

°1 °2 
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- J R(r 2 |o 6 , o 7 )i?(r 3 |o 8 , o 9 )P(o 6 , o 7 , o 8 , o 9 ), 



(21) 
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and by inserting, in the right hand side of it, the following approximation for 
the two-site probability 

P(a 1 ,a 2 )=P(a 1 )P(a 2 ). (25) 

The equation becomes then a closed equation which can then be solved by 
repeated iterations. 

Better approximations can be set up by generalizing this procedure. An 
approximation of order n is obtained as follows. From equation (|23|) we write 
down an equation for the n-site probability distribution, namely 

P(a 1 ,a 2 , — ,cr n ) = 

= E-EE w (°iK> cr' 2 )w(a 2 \a 2 , a' 3 )...w(a n \a' n , cr' n+1 )P(a[, a' 2 , ...a' n , a' n+1 ). 

(26) 

Next we use the following approximation ]TJ| \TE\ for the (n + l)-site proba- 
bility 

r>/ x P(oi,o 2 , ...,o n )P(o 2 ,o 3 , ...,o n+ i) 

P(oi,o 2 ,o 3 , ...,a n ,a n+1 ) = — - , (27) 

P{a 2 ,a 3 , ...,o„) 



and insert it in the right hand side of (|2q) . Then equation ( P6|) becomes a 
closed equation for the n-site probability distribution since 

P(o 2 , 0-3, o n ) = P(<7l, °2, °n)- (28) 
0-1 

If the order of approximation n > 4, then P(cri, <J 2 , o 3 , 04) is obtained 
from P(a 2 , cr 3 , cr n ) by summing over all Oj except a±, o 2 , o 3 , and o 4 . When 
n < 4, then for n = 1, 

P(o 1; o 2 , o 3 , o 4 ) = P^)?^)?^)/'^), (29) 

for n = 2, 

P(o 1; o 2 , o 3 , o 4 ) = p^pfo) ■ ( 30 ) 



and for n — 3, 



D / \ P(o-i,a- 2 ,a 3 )P((T2, 0-3,0-4) , Q1 v 

P(Oi, 2 , O3, O4) = — r . (31) 

P(cr 2 ,0- 3 ) 
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4 Results 



The renormalization transformation defined in the previous sections may be 
sought of as a trajectory in the space spanned by the parameter pi, and 
P2- It can be easily checked that w(l|0,0) = implying that a system with 
one absorbing state renormalizes into another one with the same property. 
Within the space of parameters p\ and P2, the renormalization transformation 
may be seen as a mapping (pi,p2) — > {v'liV^)- Given u>(l|10) = io(l|01) = pi, 
w(l|ll) = P2, and w(l|0, 0) = we obtain p[ = iu(l|10) and p' 2 = w(l|ll) 
by using equation (p~9|) . It can be easily checked that w(l|0, 0) = implying 
that a system with one absorbing state indeed renormalizes into another one 
with the same property. 

We have iterated equation flTS|) numerically using approximations of order 
n = 1, 2, 3, 4, 5 and 8. For all orders of approximation that we used we have 
found, in the plane (pi,P2), two fully attractive fixed points (0, 0) and (1, 1), 
one fully repulsive fixed point (1/2, 1) and one nontrivial fixed point (p*,^)- 
The attractive fixed point (0, 0) is related to the absorbing state whereas the 
other attractive fixed point (1, 1) is related to the active state. Almost all 
trajectories are attracted to either one or the other of these two points as 
can be seen in figure 1. The base of attraction are separated by a line, which 
should be identified as the critical line of the Domany-Kinzel model. The 
separatrix hits the line P2 = 1 at p\ = 1/2, which is a fully repulsive fixed 
point, and the line P2 = at p\ = p c . 

The nontrivial fixed point is located over the separatrix line and is attrac- 
tive along a direction parallel to the separatrix and repulsive in the direc- 
tion perpendicular (a hyperbolic fixed point). According to the RG theory 
the critical exponent u±, related to the correlation length, is obtained from 
the eigenvalue A along the repulsive direction through u± = In b/ In A. The 
eigenvalue A is the largest eigenvalue of the Hessian matrix related to the 
transformation (^1,^2) — > {v'iiV'2) calculated at the nontrivial fixed point. 
The location of the critical line (the separatrix) is in fair agreement with 



previous results coming from numerical simulation [13|, as can be seen in 



figure 2. The critical value p c = 0.792 of p\ at p 2 = differs from the value 
p c = 0.799 ± 0.002 obtained from numerical simulations in about 1%. 

Table 1 shows the values of the critical point p c (along P2 = 0) for several 
values of the order of approximation up to n = 8. For n > 4 the figures are 
very close to the value obtained from simulations. In the same table it is 
shown the values of the nontrivial fixed point together with the eigenvalue A 
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and the critical exponent v±_ = In 2/ In A. All figures for the eigenvalue A were 
obtained numerically since it is very difficult to calculated the derivatives of 
the Hessian matrix related to the transformation (pi,P2) — * (jp'itP'q)- The best 
value of the critical exponents, that were obtained from the approximation of 
order n = 8, is u± = 1.013 which should be compared to u± = 1.100 ± 0.005 
obtained from transfer matrix methods J^. 

5 Conclusion 

We used a real space RG to obtain the critical behavior of the Domany-Kinzel 
cellular automaton. The phase diagram is the same as the one obtained 
from simulation. The critical line, with the exception of the point (1/2, 1), 
is dominated by just one nontrivial fixed point implying, within the ideas of 
RG, that all critical points along the line are in the same universality class, as 
expected. The location of the critical line is in good agreement with previous 
results, the value of p\ at p 2 = differing in about 1% from the value obtained 
from numerical simulations. The best value of the critical exponents is in 
fair agreement with the values obtained from transfer matrix methods and 
from numerical simulations. A way of improving these figures is to use a 
renormalization scheme for which the renormalization factor b is less than 
two. This could be accomplished, for instance, by renormalizing a block of 
three sites into a block of two in which case we would have b = 3/2. Another 
way of getting better values is to use a RG space with more parameters 
instead of just two as we have done here. 
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Table caption 



Table 1. For each value n of the order of approximation it is shown 
the critical point p c (along p 2 = 0), the nontrivial fixed point (p*,^), the 
eigenvalue A corresponding to this fixed point along the relevant direction 
and the critical exponent u± = In 2/ In A, related to the correlation length. 
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Table 1 



n p c p\ 

1 0.745782 0.595189 

2 0.778384 0.598983 

3 0.788641 0.591302 

4 0.792203 0.591732 

5 0.791504 0.591872 
8 0.791654 0.591930 



P* 2 A u ± 

0.826419 2.118 0.924 

0.927874 2.144 0.909 

0.949335 2.034 0.976 

0.946004 2.000 1.000 

0.945275 1.988 1.009 

0.944997 1.982 1.013 
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Figure Caption 



Figure 1. Renormalization group trajectories in the parameter space 
{P11P2) °f the Domany-Kinzel cellular automaton, obtained by using the 
approximation of order n — 4. Almost all trajectories flow either to one of 
two fully atractive fixed points. One of them is related to the active state 
(open circle) and the other to the absorbing state (origin of coordinates). 
The two base of atraction are separated by a line, identified with the critical 
line of the model, where the hyperbolic nontrivial fixed point (full circle) and 
a fully repulsive point (star) are located. 



Figure 2. Phase diagram of the Domany-Kinzel probabilistic cellular 
automaton showing the frozen absorbing state and the active state. The 
squares are results coming from numerical simulations [13] and the line is 



the result obtained from the present real space renormalization group using 
the approximation of order n = 4. 
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Figure 2. Tome and de Oliveira 
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